Finite temperature stability of a trapped dipolar Bose gas 
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We calculate the stability diagram for a trapped normal Bose gas with dipole-dipole interactions. Our study 
characterizes the roles of trap geometry, temperature, and short-ranged interactions on the stability. We predict a 
robust double instability feature in oblate trapping geometries arising from the interplay of thermal gas saturation 
and the anisotropy of the interaction. Our results are relevant to current experiments with polar molecules and 
will be useful in developing strategies to obtain a polar molecule Bose-Einstein condensate. 
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Introduction: There has been phenomenal recent progress 
in experimental efforts to produce quantum degenerate gases 
with dipole-dipole interactions (DDIs) QtS]. The interest in 
these systems is being driven by a broad range of proposed 
applications from condensed matter physics to quantum infor- 
mation, e.g. see fl|5l]. For bosonic dipolar gases the majority 
of theoretical work has been undertaken for the T = case 
described by the Gross-Pitaevskii equation (e.g. see [0 | 
motivated by beautiful experiments with 52 Cr condensates | 
On the horizon is the realization of samples of degenerate po- 
lar molecules which have DDIs several orders of magnitude 
larger than those of the atomic systems. Currently bosonic 
Jit] and fermionic KRb molecules have been produced in 
their rovibrational ground state and effort is now focused on 
cooling these to quantum degeneracy. 

In addition to being a long ranged interaction, the DDI is 
also anisotropic with an attractive component. Thus an im- 
portant consideration is under what conditions the system is 
mechanically stable from collapse to a high density state. The- 
oretical studies have been performed for the (T — 0) con- 
densate, in which this stabilization arises from the quantum 
pressure of confinement (also from repulsive short range in- 
teractions ifioll ). Good agreement has been obtained between 
such theories and experiments with 52 Cr lUlll . The stability of 
dipolar condensates is considerably richer than that predicted 
for a gas with attractive contact interactions. Notably, due to 
the DDI anisotropy stability depends strongly on the geometry 
of the trapping potential ||12ill3ll . We note that detailed studies 
reveal rich nonlinear behavior (e.g. spatially oscillating con- 
densate density) in parameter regimes near the boundary of 
instability related to the emergence of roton-like excitations 
11411 . An experimental study of condensate collapse dy- 
namics has also been reported IU5II . 

Stability at finite temperature, pertinent to current experi- 
mental work with polar molecules, remains much less clear. 
In particular, while there has been some work on stability of a 
normal dipolar Fermi gas lfl6t . the finite temperature bosonic 
system remains largely unexplored. In this paper we develop 
a theory for the stability of a trapped dipolar Bose gas at tem- 
peratures above the critical temperature. Our work is based 
on self-consistent meanfield calculations in which we iden- 
tify the stability regime using the density response function, 
closely related to the system compressibility. This allows us 



to quantify the roles of trap geometry, temperature, and short 
range interactions. 

Formalism: We consider a gas of bosons that interact by 
both a long-range DDI and a contact interaction characterized 
by {7int(x) = gS(x) + £/dd(x) Q3], where g = 47rafr 2 /m, 
with a the s-wave scattering length. We take the dipoles to be 
polarized along z giving a dipole interaction of 



C/dd(x) 



Cdd 1 — 3 cos 2 

4.7T Ixl 3 



(1) 



with 9 the angle between x and the z axis. The constant 
Cdd, given by Mo Mm f° r magnetic dipoles of strength /i m and 
cP/eq for electric dipoles of strength d. The atoms are con- 
fined within a cylindrically symmetric harmonic trap U(x) = 
?j [ujp(x 2 + y 2 ) + uj 2 z 2 \ , with aspect ratio A = Ld z /Lu p . 

Our primary concern is the properties of this system above 
the critical temperature for which we use self-consistent 
meanfield theory Jl7ll giving the density as 



(2) 



where m is the chemical potential, j3 = 1/fcgT is the inverse 
temperature, Ca( z ) = iZjLi ^ / 3° i s tne Bose/Fermi 
function, and AdB = h / ^/Tjrmk^r . The effective potential, 

Kfr(x) = <7(x) + 2 ff n(x) + $ D (x), (3) 

includes direct and exchange terms for the contact interaction, 
but only the direct term for the dipole interaction, given by 



* D (x) = / ds! C/ dd (x - x')n(x'). 



(4) 



The neglect of dipole exchan ge i s consistent with other work 
on finite temperature bosons IU4I1 and zero temperature stud- 
ies of fermion stability III 811 . Dipolar exchange has recently 
been included in equilibrium calculations for the fermionic 
system |[l6-i and found to be less significant than direct inter- 
actions except in near-spherical traps IU9II . More generally, 
the experience from fermion studies suggests that exchange 
effects will give rise to shifts in the stability boundaries, but 
not change the overall qualitative behavior Ilia Il8ll . We also 
neglect collisional loss, such as exothermic bimolecular reac- 
tions. This is an issue for reactive molecules such as KRb 112011 
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but is expected to be unimportant for other alkali-metal dimers 
1 2 111 being pursued in experiments, such as RbCs [22]. 

The cylindrical symmetry of the system allows an accurate 
numerical treatment using Fourier-Bessel techniques 1123 1. 
The main challenge in solving the mean field theory is that 
the numerical grids must have sufficient range and point den- 
sity to accurately represent the DDI. In practice the grids be- 
come quite large for anisotropic traps. Another requirement is 
that Eqs. ©-(IUi need to be solved self-consistently, which we 
implement using a fixed point iteration scheme. Near insta- 
bility regions, the convergence rate of this scheme generally 
decreases significantly. 

To calculate the stability phase diagram we determine the 
parameter regime where self-consistent meanfield solutions 
are obtained. In practice we see a number of signatures of 
instability of our solutions, such as the divergence of the den- 
sity (i.e. density spike) and strong dependence on the nu- 
merical grid. While such failures of convergence have been 
widely used to identify meanfield instability, for the results 
we present here we use an unambiguous condition in terms of 
the density response function diverging [24]. This divergence 
is related to the mode-softening used in calculations by Ronen 
et al. to identify dipolar condensate instability Jgt]. 

The density response function of the system, in the random 
phase approximation (RPA), is given by 

Xo(k) 



X(k) 



l + [2.g + C/dd(k)]xo(k)' 



(5) 



Here Xo(k) is the bare response function and £/dd(k) = 
Cdd(cos 2 9k — 1/3) is the Fourier transform of C/dd(x). Within 
the semiclassical approximation used in our meanfield the- 
ory the response functions also depend on position, however 
as instability occurs at trap center we take these to be eval- 
uated at x = 0. We note that Eq. (0 includes direct and 
exchange contact interactions 12511 . but only the direct dipo- 
lar interaction (also see II26|P . For stability considerations we 
take k — > along the direction which f/dd(k) is most attrac- 
tive, i.e. 9k = 7r/2. With this limit the stability region is de- 
termined by the condition 



l + (2.g-C dd /3) X o(0) >0, 



(6) 



where X o(0) = /S^ 2 ( e ^ _K * (0)1 )AdB' (e-g- see 0271]). It is 
worth noting that the k — > limit of the density response func- 
tion is proportional to the compressibility, i.e. k , = ?i~ 2 x(k— >■ 
0). Thus instability is signaled by diverging compressibility 
and hence density fluctuations of the system. 

Results: We present results for the stability regions of a 
purely dipolar gas (g = 0) in Fig. Q2 a )- We observe that the 
stability region grows with increasing A. The strong geometry 
dependence of these results arises from the anisotropy of the 
dipole interaction: In oblate geometries (A > 1) the dipoles 
are predominantly side-by-side and interact repulsively (sta- 
bilizing), whereas in prolate geometries (A < 1) the attrac- 
tive (destabilizing) head-to-tail interaction of the dipoles dom- 
inates. Geometry dependence has also been observed in the 
stability of a (T = 0) dipolar Bose condensate @|. 
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FIG. 1. (Color) Stability of a purely dipolar Bose gas. (a) Stability 
boundaries for different trap geometries, where the shaded regions to 
the right of the boundaries are stable. Filled symbols indicate self- 
consistent calculations of the stability boundary according to Eq. l[6j. 
(b) Comparison of Bose (solid lines) and Fermi (dashed lines) gas 
stability, (c) Density response functions for A = 10 along the fixed 
temperature path A indicated in (a). The classical limit of the bare 
response function, /3n(0) is shown, where n(0) is the density at x = 
0. To make our calculations independ ent of N we scale T by the 
ideal gas critical temperature T° = y/ N/ £(3)/kj /fes and use the 
interaction parameter Co = fojaj^,/ \/N, where C( a ) = Ca~(l)> 
yjh/rrujj and uj 
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A prominent feature in Fig. [TJa) is that all the stability 
boundaries terminate at the critical point with Cdd = 0. This 
can be understood because Xo(0) of a saturated Bose gas di- 
verges and the system is unstable to the attractive component 
of the dipole interaction [see Eq. ©]. The harmonic trap pro- 
vides a long wavelength cut-off that limits the divergence of 
Xo(0) ll27ll and (beyond the semiclassical approximation) will 
allow systems with small Cdd to be stable below the critical 
point. At temperatures well-above condensation thermal pres- 
sure dominates and the critical dipole strength for stability in- 
creases with temperature. 

It is interesting to contrast the behavior to that of an 
equivalent system with Fermi statistics, for which n(x) = 
A dBC 3 ~/ 2 (e^^~ VaH W]). The stability regions, identified us- 
ing an analogous procedure to the Bose gas, are shown in 
Fig. |TJb). At high temperatures both systems exhibit simi- 
lar stability properties. The systems are distinctly different at 
low temperatures (note T F w 1.93T") with the Fermi system 
being stabilized by degeneracy pressure. We note that Fermi 
stability calculations using the same theory have been carried 
out in ifish f or T = 0. 

A striking feature of the oblate system in Fig. [TJa) is that 
the stability boundary bends back on itself so that the sys- 
tem is only stable for moderate values of Cdd- This feature, 
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which we refer to as double instability, first arises for mod- 
erate anisotropies of the confining potential (A > 2), but be- 
comes more prominent as A increases. The physical origin of 
the double instability can be understood by considering sys- 
tem properties along the vertical line marked A in Fig. Q3 a ), 
and by noting that in the purely dipolar case the stability con- 
dition (O reduces to CidXo(O) < 3. For the lowest values of 
Cdd the system is saturated (n ps 2.612/A^ b and \o — > °°) 
and unstable for any attractive interaction. However, since the 
average effect of the dipolar interaction across the whole cloud 
is repulsive for the oblate trap, the effect of increasing Cdd is 
to decrease the central density. Any decrease in density from 
saturation causes a rather large decrease in xo [see Fig.Qlc)] 
and the system becomes quite stable (i.e. \ takes a moderate 
value). As Cdd is further increased the system eventually be- 
comes unstable due to interactions (noting \o is rather small 
at this upper instability as the system is far from saturation) 
[Fig-Ec)]- Thus the lower instability arises from saturation 
(divergence of the bare response function) while the upper in- 
stability is driven by the large interaction strength. Because 
the normal Fermi gas cannot saturate, the double instability 
feature cannot occur for this system [c.f. Fig. 03b)] 




FIG. 2. (Color) Stability for various values of the contact interaction 
with (a) A = 1 and (b) A = 10. Shaded regions to the right of 
the boundaries are stable. Symbols indicate where self-consistent 
calculations determine the stability boundary according to Eq. (|6). 
Stable condensed regions are indicated for cases with g > (for 
g = condensation occurs at Cdd = with T c = T®, and for g < 
the system is unstable prior to condensation). Other parameters as in 

F5g.m 

In Fig.|2]we investigate the effect of contact interactions on 
the stability diagram for two different trap geometries. In both 
cases we note that for repulsive contact interactions (g > 0) 
the stability region is increased over the purely dipolar gas, 
while for the attractive case (g < 0) stability is reduced. These 



observations can be qualitatively understood from the stability 
condition ([6j, e.g. a positive value of g can offset the attractive 
component of the DDI. Indeed, if 2g > Cw/3 then there is no 
attractive component to the overall interaction and the system 
is stable. We see this in Fig. [2ja) and (b) where the boundary 
lines terminate at the critical temperature T c with the dipole 
interaction strength C^ d = 6g. In Fig. [2] we schematically 
indicate the stabilit y re gions for the condensed phase below 
T c using this result [28]. 

Finally we comment on the effect of contact interactions 
on the double instability feature that occurs in the oblate trap 
[Fig. |2jb)]. Attractive contact interactions make the feature 
more prominent (noting this case cannot stably condense). 
The double instability region gets smaller for moderate val- 
ues of repulsive contact interactions and admits a stable con- 
densate phase. For sufficiently large values of g the double 
instability region disappears. In the A = 10 trap this occurs at 
g « 1.5C [Fig.Eb)]. 

Discussion: Our results are most significantly applicable 
to current experiments with polar molecular gases, which are 
now approaching degeneracy. To put our results in context we 
now discuss the typical interaction parameters accessible in 
the lab. While the dipole strength is dependent on the electric 
field applied, using the maximal value for KRb of d ps 0.57 
D we have C dd ~ 4.76 x 4ttC Q. For RbCs the dipole 
moment is expected to be about twice that of KRb 1 30] . Much 
less is understood about the contact interactions of the molec- 
ular systems, although some progress on understanding the 
s-wave properties has been made in Ref. ll20ll . 

In contrast, atomic systems interact with typically much 
weaker magnetic dipoles. For comparison (using the same 
trap and particle number) the parameters of 52 Cr would be 
Cdd = 0.0117 x 47rC and g = 0.325C (which can be var- 
ied using a Feshbach resonance 111 111 ). This value of Cdd is 
rather small on the scale of our phase diagrams [Figs. [Tja) 
and|2a)-(b)], and thus instability above T c is not a concern 
for this system even in the absence of any contact interac- 
tion. Furthermore, since g ^> Cid/6 this system is stable 
below T c . As noted earlier, for the trapped system the bare 
compressibility gets large but does not diverge at T c as the 
residual quantum pressure provides some stability. An in- 
vestigation of condensate stability was performed using the 
Gross-Pitaevskii equation 12]. There it was found that stability 
depends on A in a complicated manner but mostly the stabil- 
ity boundary was seen to increase with increasing A, reaching 
only Cdd ~ 0.006 x 4ttC for A = 10 and N = 10 5 . This also 
allows us to conclude that while quantum pressure is an im- 
portant consideration for the small atomic dipoles it is rather 
unimportant for the above polar molecules at their maximal 
dipole strength. 

Conclusions: In this paper we have quantified the stability 
of a normal dipolar Bose gas. Our results quantify the rich 
interplay of DDI anisotropy, trap geometry and contact inter- 
actions in determining the stability regions, and demonstrate 
the distinctive behavior of the dipolar Bose and Fermi gases. 
We have predicted a novel double instability region in oblate 
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traps, and explained how this arises from a competition be- 
tween saturation and interaction effects. 

Pivotal to our analysis has been the use of the density re- 
sponse function in the RPA, which we find to accurately pre- 
dict where our meanfield calculations become unstable and 
provide a quantitative condition for stability of the saturated 
gas. Interestingly, experimental techniques have recently been 
developed to measure the density response function (equiva- 
lently the compressibility) in a trapped atomic gas B31I1 . Using 
these techniques should furnish a broader understanding of the 
mechanisms leading to instability. 

Our predictions will be relevant to current and emerging 
experiments with polar molecules and suggest that a range of 
strategies including the use of highly oblate traps, reduction 
of the dipole strength, or increasing the contact interaction 
strength will be necessary to have a stable pathway to cool 
the system to a Bose-Einstein condensate. Additionally, our 
prediction of a novel double instability feature occurs in oblate 
traps currently favored in experiments (e.g. see (13211 ). 

There are several extensions to the theory that will be con- 
sidered in future work: (i) Beyond semiclassical treatment 
of modes to account for the quantum pressure effects below 
T c \ (ii) The inclusion of dipolar exchange interactions; (iii) 
Beyond meanfield treatment to account for correlation in the 
strong dipole regime. Classical field methods provide an av- 
enue of investigation that accounts for many of these consid- 
erations and should be applicable at and below T c 13311 . 
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